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The distribution function of envelope voltage for short samples of atmospheric radio 
noise as received by a communications receiver in the VLF range always shows a marked 
departure from that obtained for Gaussian noise. If it is considered that this departure is 
caused by strong noise pulses which do not overlap in time, the effect of changes in the 
receiver bandwidth on the observed distribution function can be deduced by consideration 
of the changes in the receiver impulse response. A transformation can be obtained which 
gives an excellent approximation to the change in a mathematical representation of the 
distribution function in the range of probabilities below 1 percent. Empirical relationships 
are suggested which give useful estimates of the change in the distribution function over 
the total range of probabilities. 



1. Introduction 

In the VLF range, atmospheric radio noise is one 
of the important factors in a system design. If all 
of the manmade interference is controlled by alloca- 
tion procedures and good engineering, the atmos- 
pheric noise provides the ultimate background dis- 
turbance from which the desired signal must be 
separated [CCIR Rpt. 65, ITU, 1957; Watt, Coon, 
Maxwell, and Plush, 1958]. Studies of the noise 
must be of a statistical nature, bub there are problems 
which differ, for instance, from those of statistical 
studies of thermal noise, because the atmospheric 
noise is a nonstationary process. Since the statistics 
of the process change with time, the accuracy of 
measuring statistical parameters cannot be in- 
definitely increased by increasing the length of time 
of the measurements. Experience has shown that 
when measuring the amplitude probability distri- 
bution of the noise envelope, samples of noise 10 to 
20 min in length are short enough to avoid difficulties 
due to the nonstationary characteristics of the 
process, but are long enough to give useful informa- 
tion for system design purposes [Hoff and Johnson, 
1952; Watt and Maxwell, 1957]. 

One of the problems which occurs in system design 
is determining the percentage of time that the 
envelope of the noise will exceed the signal level, and 
how this varies as the receiver bandwidth is changed. 
In principle this can be calculated precisely if 
enough is known about the statistics of the noise. 
An applicable procedure, as described by Widrow 
[1957], is to consider a sampled version of the input 
noise for which the joint probability distribution of 
the samples can be determined, and to use this to 
calculate the probability distribution of a sampled 
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version of the filter output. If the sampling rate at 
the input is high enough, the probability distribution 
of the output samples will accurately represent the 
probability distribution of the continuous output. 
The difficulty with this procedure is that the interval 
between samples at the input must in general be 
short relative to the correlation time of the noise 
[Ragazzini and Franklin, 1958], so that a very high 
order joint probability density function must be 
calculated for the filter output. The purpose of this 
paper is to present a technique which requires much 
less computation, but gives an accuracy commen- 
surate with the statistical knowledge of the noise. 

2. Noise Characteristics 

Figure 1 shows two samples of measurements of 
atmospheric noise envelope distributions, one of rela- 
tively high dynamic range conditions and the other 
of noise of relatively low dynamic range; the dotted 
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Figure 1. Typical measurements of atmospheric noise envelope 
distribution. 
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line represents the measurement which is obtained 
for thermal noise. At the low voltages which have 
a high probability of being exceeded, the atmospheric 
noise measurements parallel those for thermal noise; 
this is characteristic of a phenomenon consisting of 
a large number of overlapping pulses, no one of which 
contributes a significant proportion of the total 
energy. At the high voltages, however, which have 
a low probability of being exceeded, the measure- 
ments depart radically from the shape of the thermal 
noise curve, in the form of a much higher probability 
of exceeding the voltage levels in the range of prob- 
abilities of 0.3 and below. Oscillograms of atmos- 
pheric noise show that the high voltage levels are 
attained by occasional noise pulses which greatly 
exceed the general level of the noise [Yuhara, Ishida, 
and Higashimura, 1956] and that the pulse shape is 
that of the impulse response of the receiver. These 
strong pulses occur infrequently and two such pulses 
very rarely overlap. 

3. Probabilities Associated With Impulses 

The effect of bandwidth changes on the probabili- 
ties associated with noise impulses can be calculated 
by considering the effect of bandwidth changes on 
the usual impulse response of the filter. The effect 
of increasing the bandwidth is to cause the impulse 
response to become shorter in time, and higher in 
amplitude [Guillemin, 1956]. More specifically, if 
the network characteristics are changed in such a 
way that the shape of the impulse response is pre- 
served while the bandwidth is multiplied by a factor 
B, the voltage scale of the impulse response will be 
multiplied by a factor B and the time scale by a factor 
1/B, as illustrated in figure 2 for B = 2. This rela- 
tionship also holds for the envelope of the response 
of a band-pass circuit whenever the circuit possesses 
a low pass analog, as is normally the case with 
receiver circuits [Aigrain, Teare, and Williams, 1949; 
Guillemin, 1953]. Further, the probabilities asso- 
ciated with the impulse response are determined by 
the time scale; for example, the probability that a 
voltage level of one will be exceeded by the narrow 
band response is determined by the ratio of the time 
r to the total observation time. Because of the 
relationship between changes in the voltage and time 
scales, the doubled bandwidth impulse response 
exceeds a voltage level of two for a time of r/2, 
which represents precisely one-half of the probability 
corresponding to a time r. 

Thus the amplitude distribution of a short sample 
of noise which contained only one impulse would be 
determined by the time and voltage scales of that 
impulse response; if a particular voltage v is measured 
to have a probability P of being exceeded when a 
single impulse is observed through a network of band- 
width Aco!, and simultaneously the impulse is ob- 
served through a network of bandwidth Aoo 2 = BAu 1 , 
the voltage Bv will be measured to have a probability 
(l/B) P of being exceeded. This relationship also 
holds for a noise sample which consists of a number 
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Figure 2. Envelope of band-pass impulse response for two 
banalwidths. 



of impulses, provided that the response of the net- 
work to one impulse does not overlap the response 
to another. 

These requirements are never exactly met in prac- 
tice because of the exponential decay of physical 
networks and because of the presence of thoroughly 
overlapping low level noise. However, the high 
voltage end of the distribution of VLF atmospheric 
noise appears to represent an excellent approximation 
to the requirements. Figure 3 shows the high volt- 
age portion of two representative sets of measure- 
ments at two bandwidths and also shows the results 
of translating the narrow-band measurements by the 
ratio of bandwidths to expected measurements at 
the wider bandwidth. It is evident that for both the 
low and high dynamic range noise conditions the 
points calculated from the narrow-band measure- 
ments agree with the measurements at the wider 
bandwidth for probabilities less than 0.01. 

4. Mathematical Representation of the Noise 

In dealing with a non-Gaussian random variable, 
it has proved profitable for some purposes to consider 
it as being generated as the output of a nonlinear 
resistance network which has a Gaussian input, even 
though the actual phenomenon may be very much 
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Figure 3. Atmospheric noise distributions showing effect of translating measurements made at a narrow band- 
width to expected values at a wider bandwidth. 



different from this [Smith, 1959a and b]. In the 
situation being considered here where the envelope 
of the noise is the quantity of interest, we can 
represent the envelope of the actual noise as being 
obtained from an envelope detector having a Gaus- 
sian input, followed by a nonlinear resistance net- 
work. The justification for doing this is simply that 
the resulting probability distribution can represent 
the probability distribution of the actual phenomenon 
quite accurately, and the model proves to be a con- 
venient one for the type of mathematical manipula- 
tions which are required. 

A nonlinear network which can produce the correct 
probability distribution has a voltage transfer char- 
acteristic which is a three term polynomial of the 
form 



v=a x y+a 2 y 2 +a z y b 



(1) 



where a h a 2 , a 3) and b are chosen to match the 
atmospheric noise distribution under consideration; 
the probability that a voltage level v is exceeded at 
the output of the network is given by 



P=e-y\ 



(2) 



The variable y may conveniently be considered just 
as a parametric variable, although it is the envelope 
voltage of a Gaussian variable having a variance of 
one-half. 



Each term of the polynomial, if considered sepa- 
rately, plots as a straight line of log-log of proba- 
bility versus logarithm of voltage, which are the 
coordinate scales used on the graphs of distributions 
in figures 1 and 3. The first term represents a dis- 
tribution having complete overlapping, such as 
thermal noise; the second and third terms represent 
departures from this distribution. The third term 
of the polynomial, a 3 y b , is the dominant term at the 
high voltage, low probability end of the curve, when 
the phenomena forming the distribution appear to 
be strong, nonoverlapping pulses. 

5. The Effect of Bandwidth Changes on the 
Mathematical Model 

Since the third term of the polynomial is consid- 
ered to represent nonoverlapping impulses, the effect 
of bandwidth changes can be obtained by the proce- 
dure described in section 3. That is, for an increase 
in bandwidth by a factor B, any point on the original 
line at voltage v and probability P is transformed 
to a voltage Bv at a probability P/B. After this is 
done, it is found that the points no longer lie on a 
straight line, but along a line having a slight curva- 
ture, as illustrated in figure 4 for an exponent of 8 
and a bandwidth increase of ten times. The curva- 
ture introduced in this way is so slight, however, 
that the curve may be replaced by a line tangent to 
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Figure 4. Curvature introduced by bandwidth transformation. 
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Figure 5. Tangent approximation to curve. 
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it at the point where P=e~ 8 ; as is apparent from 
figure 5, this is an excellent approximation. 

It is possible to derive a general formula for the 
effects of bandwidth changes on the exponent and 
the coefficient of the third term of the polynomial. 
For a bandwidth increase by a factor of B, the expo- 
nent b w for the wider bandwidth is obtained from 
the exponent b n for the narrow bandwidth by the 
relationship 



&«,= 



S-L n B 



ts b n 



(3) 



where L n B is the natural logarithm of B. The co- 
efficient is obtained from 



&3w — QZnB 



[(8-L n B)$(*-^ B }~]z (4) 



where a 3w and a Zn are the coefficients for the wide and 
narrow bandwidths respectively. The derivation of 
these relationships is given in the appendix. 

For the other terms of the probability function the 
effects of bandwidth changes are not as apparent as 
for the third term which describes the behavior at 
high voltages. The second term of the polynomial 
is important in the moderate voltage range where 
some overlapping of noise pulses occurs. For the 
values of exponents and coefficients which are re- 
quired to match the mathematical model to the 
observed atmospheric noise conditions, it happens 
that the second moment of the distribution is almost 
entirely determined by the second and third terms 
of the polynomial. Since for normal communica- 
tions bandwidths the total noise power received is 
directly proportional to the bandwidth, it is possible 
to adjust the second term to obtain the correct 
behavior for the second moment of the distribution. 



The procedure is to calculate the changes in the a 3 
coefficient and exponent b as described above, and 
to insert these into the expression for the mean 
squared voltage of this form of distribution, which is 

{v*)=a\+a\T ( b -^)+alT(b+l)+2a 1 a 2 T (^±Z) 

+ 2a 1 a 3 r (^)+2a 2 a 3 r (^) (5) 

where (v 2 ) is the mean squared voltage. 

For a bandwidth increase by a factor of B, the mean 
squared voltage for the wide bandwidth is B times 
the mean squared voltage at the narrow bandwidth. 
With this and a value for a i} equation (5) can then be 
solved for the coefficient a 2 . 

The first term of the polynomial, involving a u 
represents a situation where complete overlapping 
occurs. The slope of a distribution represented by 
only this term does not change with bandwidth, but 
only its amplitude, as long as the complete over- 
lapping is maintained. If this were the only term 
in the distribution, its amplitude would change 
precisely as the square root of the bandwidth. 
Since there are other terms present, however, there 
will be an interchange of energy between the various 
portions of the distribution, and the net effect will 
be that the amplitude of the first term will change 
more slowly than this; an increase according to the 
four-tenths power of the bandwidth appears to fit 
the available experimental data quite well. 

These procedures work quite well with noise 
samples of low and moderate dynamic range, but 
an error appears when used with noise of high dy- 
namic range, that is, values of b of about 11 or 
greater. This occurs because the mathematical 
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model being used is not sufficiently complex. With 
the very Sigh dynamic range conditions, a sub- 
stantial contribution to the integral for the mean 
squared voltage occurs at probabilities below 10" 6 ; 
in this region the tangent line and the curve repre- 
senting the distribution diverge significantly, causing 
considerable error in the evaluation of the mean 
squared voltage. 

Under the high dynamic range conditions it seems 
appropriate to make the same transformation on the 
second term as on the third. Modification of the 
coefficient by some factor to compensate for the 
partial overlapping might be considered, but also 
considering the nature of the experimental data, 
that is, short samples of a nonstationary process, 
and also the relative unimportance of the second 
term, any refinement in the transformation seems 
unjustified. 

6. Discussion 

These somewhat arbitrary operations on the first 
and second terms, together with the accurate trans- 
formation for the high voltage end of the distribution 
where the third term is dominant, provide a method 
of describing what changes would be observed in 
the polynomial representation of the distribution of 
atmospheric radio noise amplitude when observed 
through different bandwidths. The results of carry- 
ing out the complete transformation on representa- 
tive samples of high and low dynamic range noise 
are shown in figure 6 and figure 7. For the low 
dynamic range noise in figure 6, the a 2 coefficient 
was determined from the rms voltage relationship; 
for the high dynamic range noise in figure 7, it was 
determined by the transformation procedure for 
nonoverlapping pulses. It is apparent that those 
procedures provide good estimates of the distri- 
butions, ami are useful and practical tools for use 
when more rigorous methods are not justified. 
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Figure 6. Low dynamic range atmospheric noise distribution 
measured at a wide bandwidth, and estimated distribution 
from simultaneous narrow bandwidth data. 
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Figure 7. High dynamic range atmospheric noise distribution 
measured at a wide bandwidth, and estimated distribution 
from simultaneous narrow bandwidth data, 
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7. Appendix: Variation of High Voltage Co- 
efficient and Exponent 

Consider a distribution having a form expressed by 



Vi=ay v 



Pi = e~ 



(6) 



when it is observed through a particular bandwidth. 
If this distribution is composed of nonoverlapping 
impulses, then when observed through a bandwidth 
B times larger, the voltages and probabilities will be 
related by the ratio of bandwidths according to 



pJPi 

Vi=Bvi 



(7) 



where P x and e?j are obtained at the first bandwidth 
and P 2 and v 2 at the second. From (6) and (7) we 
can write 

]og(-£nP,)-flog(£) 

log(-L ft BP 2 )=?log(^) 



(8) 



We wish to develop an approximation to the rela- 
tionship between v 2 and P 2 in the region of high volt- 
ages and low probabilities which will have the form 



v,=aBC 1 y K * 

P 2 =e-** 



(9) 
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or 



log(-ZnP 2 ) = Al og (j0 



(10) 



where C ] and C 2 are constants to be determined. 
They may be evaluated by selecting two points on 
the P 2 versus v 2 relationship where the approxima- 
tion is to be exact. Choosing the points P 2a , v 2a , and 
P 2 p } v 2 p, we obtain from (8) 



v 2a =aB(—LnBP 2a ) 2 

± 
v a fi=aB(-LnBP 2fi )*. (11) 

Now in the approximation, 

log(-LnP 2a ) = Ai og (j^_J 

iog(-i^)=A log (_^.) (12) 

(7imay be eliminated from equations (12), obtaining 



and from (11), 



V2a/ LnBP 2a \* 
v 2 p \LnBP 2 8j 



o 7 

8/ 



(14) 



and from (13) and (14) we can obtain 
. ( LnBP 2a \ 

" , (LnP 2o \ ' 



C 2 - 



(15) 



For the evaluation of C u from (11) and (12) we can 
obtain 



Cr- 



{-LnBP 2a )* _ {-LnBP 2B y 

bC 2 hc£ 

(-LnP 2a y (-LnP*)~ r 



(16) 



Going back to the polynomial term itself, we can 
change the third term to read a$ w y hw , where 



and 



b w =b n C 2 

Uzw ==a ZnBCi 



(17) 



dzn and b n being the values observed at the narrow 
bandwidth. 

A limiting value of C 2 is obtained as P 2a approaches 
P 2/ 3 and this represents the approximation where 
the straight line is tangent to the curve. This 
limiting value may be found by application of 



THospitaFs rule, and is 



c 2 = 



LnP, 



23 



LnBP# 



(18) 



If we use the tangent approximation and evaluate for 
p 2a= p 2/3==e -8 ==0 .000335 



we obtain 



and 



b w =C 2 b n = 



^3; 



S-LnB " n 



n f 8-LnB \ 
,-<hn* ^ 8 c 2 J 



(19) 



(20) 



Figure 5 shows this tangent line for an assumed b n of 
8 and B of 10, and also shows some points from the 
curve in figure 4. The difference between the 
straight line and the points from the curve is zero, 
of course, at P= 0.000335, and increases slowly to 
about 1 db at P=0.005 and P=10" 6 ; it is evident 
that there is no significant difference between using 
the straight line and the curve, and the procedure 
using the tangent line provides a simple and accurate 
estimate of the changes in the high voltage end of 
the distribution curve. 
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